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We study the effect of localisation on the propagation of a pulse through a multi-mode disordered 
waveguide. The correlator {u{uji)u* {uj2)} of the transmitted wave amplitude u at two frequencies 
differing by Slo has for large Slo the stretched exponential tail oc exp(— -y/ tbSlj/2). The time constant 
td ~ L^/D is given by the diffusion coefficient D, even if the length L of the waveguide is much 
greater than the localisation length ^. Localisation has the effect of multiplying the correlator by a 
frequency- independent factor exp(— which disappears upon breaking time-reversal symmetry. 

PACS numbers: 42.25.Dd, 42.25.Bs, 72.15.Rn, 91.30.-f 



The frequency spectrum of waves propagating through 
a random medium contains dynamical information of in- 
terest in optics acoustics Q, and seismology A 
fundamental issue is how the phenomenon of wave locali- 
sation Q affects the dynamics. The basic quantity is the 
correlation of the wave amplitude at two frequencies dif- 
fering hy 6ll}. a recent microwave experiment by Genack 
et al. measured this correlation for a pulse transmitted 
through a waveguide with randomly positioned scatter- 
ers. The waves were not localised in that experiment, 
because the length L of the waveguide was less than the 
localisation length ^, so the correlator could be computed 
from the perturbation theory for diffusive dynamics [^. 
The characteristic time scale in that regime is the time 
Td = L'^ / D it takes to diffuse (with diffusion coefficient 
D) from one end of the waveguide to the other. Accord- 
ing to diffusion theory, for large 5uj the correlator decays 
oc exp(— y/t£)(5(jj/2) with time constant td- 

What happens to the high-frequency decay of the cor- 
relator if the waveguide becomes longer than the local- 
isation length? That is the question addressed in this 
paper. Our prediction is that, although the correlator is 
suppressed by a factor exp(— the time scale for the 
decay remains the diffusion time td — even if diffusion 
is only possible on length scales ^ L. The exponential 
suppression factor disappears if time-reversal symmetry 
is broken (by some magneto-optical effect). Our analyti- 
cal results are based on the formal equivalence between a 
frequency shift and an imaginary absorption rate, and are 
supported by a numerical solution of the wave equation. 

We consider the propagation of a pulse through a dis- 
ordered waveguide of length L. In the frequency do- 
main the transmission coefficient tnm{^) gives the ra- 
tio of the transmitted amplitude in mode n to the in- 
cident amplitude in mode m. (The modes are normal- 
ized to carry the same flux.) We seek the correlator 
C{6uj) = {tnm{uj + 5u))tl^{uj)) . (The brackets (• • •) de- 
note an average over the disorder.) We assume that 
the (positive) frequency increment 5uj is sufficiently small 



compared to lo that the mean free path / and the number 
of modes N in the waveguide do not vary appreciably, 
and may be evaluated at the mean frequency uj 0|. We 
also assume that I ^ c/lo (with c the wave velocity). The 
localisation length is then given by ||] ^ = {(3N -f 2 — [3)1, 
with [3 — 1{2) in the presence (absence) of time- reversal 
symmetry. For 3> 1 the localisation length is much 
greater than the mean free path, so that the motion on 
length scales below ^ is diffusive (with diffusion coeffi- 
cient D). 

Our approach is to map the dynamic problem with- 
out absorption onto a static problem with absorption . 
The mapping is based on the analyticity of the trans- 
mission amplitude tnm{^ + ij/), at complex frequency 
Lu + iy with y > 0, and on the symmetry relation 
tnmi^+iy) ^ inmi^^^^y)- The product of transmission 
amplitudes i„m(a; -I- z)tnm{—^ + z) is therefore an ana- 
lytic function of z in the upper half of the complex plane. 
If we take z real, equal to \5u), we obtain the product 

of transmission amplitudes tnmi^^ + 5'^^)Cm('-^ ~ l^"-^) 
considered above (the difference with t„m(w -I- (5a;)i*^(w) 
being statistically irrelevant for 5u! <C w). If we take 
z imaginary, equal to i/2Ta, we obtain the transmission 
probability T — \tnm{^ + i/2''"a)P at frequency u and 
absorption rate 1/ra. We conclude that the correlator 
C can be obtained from the ensemble average of T by 
analytic continuation to imaginary absorption rate: 



CiSij) = (T) for 1/ra 



-iSuj. 



(1) 



Two remarks on this mapping: 1. The effect of ab- 
sorption (with rate 1/t*) on C{Slo) can be included by 
the substitution 1/ra — —iSuj + 1/t* . This is of impor- 
tance for comparison with experiments, but here we will 
for simplicity ignore this effect. 2. Higher moments of 
the product C = tnm{i-^ + 5<^<^)Cm('^ ~ 5'^'^) related 
to higher moments of T by (C^) = (TP) for 1/ra —iSuj. 
This is not sufficient to determine the entire probability 
distribution P{C), because moments of the form {C^C*'^) 
can not be obtained by analytic continuation [|o|. 
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To check the validity of this approach and to demon- 
strate how effective it is we consider briefly the case 
= 1. A disordered single- mode waveguide is equiv- 
alent to a geometry of parallel layers with random vari- 
ations in composition and thickness. Such a randomly 
stratified medium is studied in seismology as a model for 
the subsurface of the Earth . The correlator of the re- 
flection amplitudes K{Suj) — {r{uo + 5uo)r*{Lo)) has been 
computed in that context by White et al. lO] (in the 
limit L oo). Their result was 



K{5uj) = {2l/c)5uj / dx eyi^[~x{2l/c)6uj] 



(2) 



The distribution of the reflection probability R = |rp 
through an absorbing single-mode waveguide had been 
studied many years earlier as a problem in radio- 
engineering , with the result 

(R) - (//cTa) / dz exp[-(z - l)(//cTa)]— -. (3) 
Ji z + 1 

One readily verifies that Eqs. (Q) and (||) are identical 
under the substitution of 1/ra by —iSu. 

In a similar way one can obtain the correlator of 
the transmission amplitudes by analytic continuation to 
imaginary absorption rate of the mean transmission prob- 
ability through an absorbing waveguide. The absorbing 
problem for = 1 was solved by Freilikher, Pustilnik, 
and Yurkevich That solution will not be consid- 

ered further here, since our interest is in the multi-mode 
regime, relevant for the microwave experiments The 
transmission probability in an absorbing waveguide with 
iV > 1 is given by Q 



(T) 



iV^aSinh(L/^a) 



cxp 



3/3,1 



—) 
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(4) 



for absorption lengths = VDt^ h^ the range Z ^ ^a <C 
^. The length L of the waveguide should be ^ I, but the 
relative magnitude of L and ^ is arbitrary. Substitution 
of 1/Ta by —iSuj gives the correlator 



Ci6u;) 



l^/—iTo5u) 



NL sinh ^^itdSuo 



exp 



L 



2Ni) 



(5) 



where td = 1? j D is the diffusion time. The range of va- 
lidity of Eq. is L/^ ^ y/ToSuj <C L/l, or equivalently 
D/S^^ <^ Slli c/l. In the diffusive regime, for L <C ^, 
the correlator (||) reduces to the known result |^ from 
perturbation theory. 
For max 

(I?/L2,£)/^2-) ^ <^ c/l the decay of the 
absolute value of the correlator is a stretched exponential. 



21 



VtdSuj exp (~\/b^^~Sl3.iJ^] 



(6) 



In the localised regime, when ^ becomes smaller than 
i, the onset of this tail is pushed to higher frequencies, 



but it retains its functional form. The weight of the tail 
is reduced by a factor exp{—L/2Nl) in the presence of 
time-reversal symmetry. There is no reduction factor if 
time-reversal symmetry is broken. 
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FIG. 1. Frequency dependence of the logarithm of the 
absolute value of the correlator C{5ijj). The data points fol- 
low from a numerical simulation for N = 5, the solid curve 
is the analytical high-frequency result (^) for A'^ ^ 1 (with 
l3 — 1). The decay of the correlator is given by the diffusive 
time constant td = L /D even if the length L of the wave- 
guide is greater than the localisation length ^ = 6 Z. The offset 
of about 0.6 between the numerical and analytical results is 
probably a finite- A'^ effect. 

To test our analytical findings we have carried out nu- 
merical simulations. The disordered medium is mod- 
eled by a two-dimensional square lattice (lattice con- 
stant a, length L, width W). The (relative) dielectric 
constant e fluctuates from site to site between 1 ± (5e. 
The multiple scattering of a scalar wave ^ (for the case 
13 — 1) is described by discretizing the Helmholtz equa- 
tion [V^ -I- (a;/c)^e]\l/ = and computing the transmis- 
sion matrix using the recursive Green function technique 
JlSf . The mean free path I is determined from the aver- 
age transmission probability (Trtt^) = A^(l -I- L/l)~^ in 
the diffusive regime The correlator C is obtained by 
averaging t„m(w-l-(5w)t*„(w) over the mode indices n, m 
and over different realisations of the disorder. We choose 
= 2 (c/a)2, 5e = 0.4, leading to I ^ 22.1 a. The width 
ly = 11a is kept fixed (corresponding to A^ = 5), while 
the length L is varied in the range 400-1600 a. These 
waveguides are well in the localized regime, L/^ ranging 
from 3-12. A large number (some 10^-10^) of realisa- 
tions were needed to average out the statistical fluctua- 
tions, and this restricted our simulations to a relatively 
small value of A^. For the same reason we had to limit 
the range of Suj in the data set with the largest L. 
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Results for the absolute value of the correlator are plot- 
ted in Fig. |l| (data points) and are compared with the 
analytical high-frequency prediction for iV 3> 1 (solid 
curve). We see from Fig. |l]that the correlators for differ- 
ent values of L/^ converge for large 5lo to a curve that lies 
somewhat above the theoretical prediction. The offset is 
about 0.6, and could be easily explained as an 0(1) un- 
certainty in the exponent in Eq. (|l|) due to the fact that 
N is not ^ 1 in the simulation. Regardless of this off- 
set, the simulation confirms both analytical predictions: 
The stretched exponential decay oc exp(— •\/rD(5tj/2) and 
the exponential suppression factor exp(— L/2^). We em- 
phasize that the time constant td — / D of the high- 
frequency decay is the diffusion time for the entire length 
L of the waveguide — even though the localisation length 
5 is up to a factor of 12 smaller than L. 

We can summarize our findings by the statement that 
the correlator of the transmission amplitudes factorises 
in the high-frequency regime: C — > /i((5w)/2(^). The fre- 
quency dependence of fi depends on the diffusive time 
through the waveguide, even if it is longer than the local- 
isation length. Localisation has no effect on /i , but only 
on /2. We can contrast this factorisation with the high- 
frequency asymptotics K fsiSw) of the correlator of 
the reflection amplitudes. In the corresponding absorb- 
ing problem the high-frequency regime corresponds to an 
absorption length smaller than the localisation length, so 
it is obvious that K becomes independent of ^ in that 
regime. The factorisation of C is less obvious. Since the 
localized regime is accessible experimentally , we be- 
lieve that an experimental test of our prediction should 
be feasible. 
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